Critical Kondo destruction in a pseudogap Anderson model: 
scaling and relaxational dynamics 
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We study the pseudogap Anderson model as a prototype system for critical Kondo destruc- 
tion. We obtain finite-temperature {T) scaling functions near its quantum-critical point, by using 
a continuous-time quantum Monte Carlo method and also considering a dynamical large-N limit. 
We are able to determine the behavior of the scaling functions in the typically difficult to access 
quantum-relaxational regime {hw < ksT), and conclude that the relaxation rates for both the spin 
and single-particle excitations are linear in temperature. We discuss the implications of these results 
for the quantum-critical phenomena in heavy fermion metals. 
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Continuous zero temperature phase transitions in 
strongly correlated electronic and atomic models have 
attracted considerable attention as a new paradigm for 
addressing the universal features of correlated quantum 
systems [l| . Quantum criticality links two nearby phases 
and determines the physical properties in a large range 
of temperature and control parameters, the quantum- 
critical region, that fans out from the quantum-critical 
point (QCP). This paradigm is especially pertinent to 
the understanding of intermetallic rare earth compounds. 
The phase diagram of these heavy fermion metals close to 
the border of antiferromagnetism features a QCP, but the 
associated quantum-critical properties arc highly unusual 
when viewed from the standard description based on 
Landau's notion of order-parameter fluctuations Es- 
pecially, inelastic neutron-scattering measurements have 
shown that the dynamical spin susceptibility in the 
quantum-critical regime features a linear-in-T spin re- 
laxation rate and satisfies a frequency over temperature 
(w/T) scaling [|. Very recently. Hall-effect measure- 
ments have indicated that the single-particle relaxation 
rate in the quantum-critical regime is also linear in T . 

These dynamical scaling and relaxational properties 
provide important clues to the nature of the heavy- 
fermion QCP. Yet, theoretically, such real-frequency be- 
havior is difficult to study. Two regimes need to be distin- 
guished: the quantum coherent (w > T) and quantum- 
relaxational {lo < T) regimes [i {h and fc_B are set to 
1). Calculation methods (such as Monte Carlo simula- 
tions) typically work in the imaginary-time domain, and 
the nonzero Matsubara frequencies (u;„) are necessarily 
in the |a;„|/T > 1 regime. Extracting the behavior at real 
frequencies requires an analytical continuation, which is 
in general a numerically ill-conditioned procedure. The 
numerical renormalization group operates on the real- 
frequency axis, but it is not reliable for the quantum- 
relaxational regime at nonzero temperatures. 

In this Letter, we address the dynamical and relax- 



ational properties of the particle-hole symmetric pseu- 
dogap Anderson model in both frequency regimes. Our 
motivations to study this model are multifold. In lo- 
cal quantum criticality for heavy-fermion metals, the 
critical destruction of the Kondo effect is local in 

space, and the resulting interacting critical modes are 
manifested in local correlators which can be studied in 
quantum-impurity problems. The pseudogap Anderson 
model is the simplest impurity problem that contains the 
physics of critical Kondo destruction; it is well known 
that varying the Kondo coupling yields a QCP dMi, 
which separates a Kondo-screened Fermi-liquid phase 
from a Kondo-destroyed local-moment phase. However, 
a proper understanding of the dynamical scaling at fi- 
nite temperatures and the associated relaxational be- 
havior is not yet available even in this simplest model. 
Furthermore, the pseudogap Andcrson/Kondo model is 
relevant in a number of realistic physical settings. It 
has been invoked in the context of non magnetic impu- 
rities in cupratc superconductors [l3|- It has also been 
shown that a judicious tuning of a double quantum-dot 
system can produce a pseudogap in the effective den- 
sity of states [31 • In disordered metals, a novel phase 
has been attributed to the occurrence of local pseudo- 
gaps near the Fermi energy at local- moment sites [l9l |. 
Finally, the pseudogap Kondo model is the appropriate 
model to describe point defects in graphene (20[. 

We study the model by using a continuous-time quan- 
tum Monte Carlo approach (CT-QMC) [2lj. We deter- 
mine the full scaling functions at real frequencies and 
finite temperatures for both the dynamical spin suscep- 
tibility and single-electron Green's function. We achieve 
this by taking advantage of insights gained from exact 
calculations at real frequencies and finite temperatures 
in a dynamical large- iV limit of the model. The results in 
the large- limit motivate us to analyze the imaginary- 
time correlators in the physical N = 2 model in a 
way that uncovers the form of a boundary conformally- 
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invariant fixed point. The latter, in turn, can readily 
be analytically-continued to real frequency at finite tem- 
peratures. We establish that both the dynamical spin 
susceptibility and single-electron Green's function dis- 
play an Lo /T-scaling and contain a linear-in-T relaxation 
rate. As a by-product, we show that the CT-QMC ap- 
proach, which is based on a high-temperature expansion, 
can reach low-enough temperatures with enough accu- 
racy to resolve quantum-critical features. 

Pseudogap Kondo model in a dynamical large-N 
limit: To set the stage for the CT-QMC study, wc start 
with the SU(N)xSU(M) Kondo model in the pres- 
ence of a pseudogap in the limit of large N and M. In 
what follows, we set ?i = fc^ = 1. The Hamiltonian is 
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Here, the spin and channel indices are cr = 1, . . . , A'^ and 
a = 1, . . . ,M, respectively. The conduction electron den- 
sity of states takes the form: 
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FIG. 1: Scaling functions for the imaginary part of (a) Green's 
function ^(w, T) and (b) susceptibility for r = 0.3 and 

K = 0.5 at the critical ^ 1.54 (with Tj^ « 0.3i3 for r = 0). 
Both functions display oj/T-scaling with scaling functions $ 
obeying ^(uj/T — > 0) — > c or for Q or x, where c 7^ is a 
constant. (c),(d) the scaling functions in imaginary time. 



with 2D being the bandwidth. That this limit has a 
nontrivial QCP can be seen through the particular form 
of the perturbative (in r) renormalization group equa- 
tion [10|. In the limit of large N and M, the renormaliza- 
tion group beta function becomes /3(?') = —j{r — j-\-Kj'^), 
with j = Jk /D and M/N^^. This establishes that 
the QCP survives the largc-N limit and can be accessed 
perturbativcly. To order r, the large-N beta function is 
identical to its N= 2 counterpart [10|,ll2| suggesting that 
the universal critical scaling properties of the N ~ 2 
QCP are preserved by taking the large-N limit. In this 
limit, the local degrees of freedom are expressed in terms 
of pseudofermions fa and a bosonic decoupling field Ba, 
where Sa^a' = flfcr' ~ ^um'Q/N, and Q is related to the 
chosen irreducible representation of SU(N) [22I, [iJl. The 
large- equations are 



Gg^iivn) = 1/ Jk -^B{ivn)] 
Gj^{iuj„) = ia;„ - A - i;/(ia;„); 



(3) 



together with a constraint G/(r 0~) = Q/N |22| . 
Here, A is a Lagrangian multiplier enforcing the con- 
straint and ^0 = — (TrCrrn {t)cI^ (0))n is the noninteract- 
ing Green's function [13|. 

By solving the large-N equations in real frequencies for 
arbitrary uj and T [23[ , the full scaling functions in both, 
the quantum coherent (oj > T) and relaxational (T > to) 
regimes are obtained. At the critical coupling Jc(?'), we 
find that all the correlators display an w/T-scaling. This 
is demonstrated in Fig. [TJa) for the local singlcparticle 
Green's function [i.e., the T-matrix, Q{ui,T), associated 



with g{T) = Gf{T)GB{T)], and in Fig. (Hb) the local 
spin susceptibility x(cj,T'), which corresponds to xi''') = 

-G}{T)Gf{-T). 

A key insight from the large- A^ result is that the scaling 
functions contain more information beyond uj/T scaling 
per se. They have the particular form associated with 
a boundary conformally-invariant fixed point, depending 
on r as a power law in ttT/ smiirrT) |25|. To see this, 
we obtain the imaginary-time dependence from the real- 
frequency results via 



$(t) = -V j 



exp(— Toj) 
exp(— /Scj) — 7/ 



Im($(w + iO+)), (4) 



for < r < /?. Here, 77 ± for bosonic/fermionic 
$. Figure [1] shows the (c) Green's function g{T,T) 
and (d) susceptibility x{'''iT) versus the combination 
ttT/ sin(7rTr). Both collapse on a single scaling curve 
in terms of 7rT/(sin(7rrT)) for all (low-enough) T. A 
power-law behavior for r — > l/(2r) is seen over about 
7 decades, and the exponents are compatible with those 
for the frequency dependence. 

Pseudogap Anderson model at N = 2: Guided by 
the large-N results, we turn to the scaling functions 
for Q{t,T) and xi^^T) of the particle-hole symmetric 
pseudogap Anderson model at A = 2; the low-energy 
properties of this model are identical to its pseudogap 
Kondo counterpart. To this end, we bring to bear the 
recently developed hybridization-expansion Monte Carlo 
method 2l|, |26| on a quantum-critical model. This CT- 
QMC approach involves a stochastic sampling of a per- 
turbation expansion in the host-impurity hybridization 
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FIG. 2: Dynamical local susceptibility x(''", T) versus tT for 
r = 0.4, and To = O.ID with [/ = and 17 = 0.3D (a). 
Finite-temperature scaling of the Binder cumulant B{U,T) 
as a function of U at various temperatures (b), error bars 
are obtained from a jackknife error analysis. From the in- 
tersection of the curves we determine the critical point to be 
Uc{r = 0.4)/L> = 0.085 ± 0.002. 



or a weak coupling expansion [21| 
are free of any finite-size effects [29| . 

The Anderson impurity model is defined by H 
J2cr H^f^ where 



The results 



k,(7 (7 

= {^dkdUka + H. c.) (5) 



with fika = c\^Cka, fida = d^^d^, Efc being the host 
dispersion, Vdk the hybridization, and ed the impurity 
level energy. We consider the particle-hole symmet- 
ric case where ed — ~\U, with U being the onsite 
Coulomb repulsion. The host-impurity coupling is spec- 
ified by the imaginary part of the hybridization function 
r(a;) = TT^^ |Vrffcp(5(aj - e^). As in Eq. we choose 
r(a;) = To 1^1'^ Q{D — The critical point exists only 
for < r <| 

Central to the CT-QMC approach adopted here 
is the expansion of the partition function Z = 
Tr{f^e-^M7„ exp[- /jf dr Hi''\T)]} in the hybridiza- 
tion term [21|. We measure the single-particle Green's 
function Gct{t) = {Trda{T)d'l), the local spin susceptibil- 
ity x{t) {TtS,{t)Sz{0)) and powers of the local mag- 
netization (Mi;-) = {{^ drS^ir))'') where S^(t) = 
i[?T--f-(T) — ^^(t)]. The static susceptibility is obtained 

from xi'-^ = 0) = lo '^'^ Thermalization can be 

traced by (ud) which obeys (ud) = 1 in the particle-hole 
symmetric model. We also performed a binning analysis 
and obtained the integrated autocorrelation time which 
increases with decreasing temperature but turned out to 
be small (compared to the number of measurements) at 
all temperatures. For the lowest temperature consid- 
ered (PD = 9, 000) we performed 800, 000 Monte Carlo 
steps for thermalization, 1,500 Monte Carlo steps be- 



tween each measurement and 18, 750 measurements. A 
Monte Carlo step consists of an attempt to remove, insert 
and shift a segment as described in Reference [2lj . 

By varying U we can tune the model through a QCP. 
Correspondingly, Figure[5l^a) shows that the large-/? limit 
of x(''' = /3/2,/3) vanishes for small U (Kondo-screened 
phase) and is equal to the Curie constant for large U 
(Kondo-destroyed local-moment phase). To accurately 
determine Uc{r) we apply finite-temperature scaling to 

the Binder cumulant, B{U,T) = where 1/T = (3 

plays the role of the system size. We find swap moves 
between up and down spin segments [sot are necessary 
to accurately measure the Binder cumulant; for the re- 
sults in Figure[5Jb) we performed a swap move every 100 
measurements. The nature of the intersection of the data 
in Figure [2jb) implies that the phase transition is con- 
tinuous, from the location of the intersection we obtain 
the critical value of U. In the quantum-critical regime 
the static local susceptibility displays an anomalous r- 
depcndent exponent; we find 



x(r, Uc 



0.4) ~ T" 



(6) 



with X ~ 0.68(3) in good agreement with numerical 
renormalization group results |12|. 

We now discuss the finite-temperature dynamical scal- 
ing properties of Q{t,T) and xi^iT). Guided by 
the large-iV results, we plot them as functions of 
(7rT)/sin(7rTT) in Figure [3l Excellent scaling collapse 
is observed over about two decades, for all temperatures 
in the scaling regime. Wc reach an important conclusion: 

Xcrit{T,T) = ^{-^, — ^) ~ {-^, — ^) , (7) 



^sin(7rrr) ' 



'sin(7rrr) ' 



for r-i < T^, Figure [3l;b). Since < 1 - a; < 1, the 
results for x(JtT) imply that the order-parameter sus- 
ceptibility shows w/T-scaling. A similar conclusion ap- 
plies to Q{t,T), as seen in Figurf|3Ka). Our results yield 
CJ(t, T — !• 0) ^ , with the exponent 5 = 1 ~ r, which 



is believed to be exact [16[. The fact that 2(5 7^ 1 — x 



signifies the importance of vertex corrections and in part 
reflects the interacting nature of the QCP (see below) . 

The boundary conformally-invariant form of x ^^'^ 
Q immediately imply that their dependence on real 
frequency satisfies lo/T scaling and that their relax- 
ation rates, defined in the quantum-relaxational regime, 
are linear in T. Expressed in terms of Tm = 
i(91nM(a;,r)/9a;|t^=o)^^ for a correlator Af, the relax- 
ation rates = aT and Tg = bT, where a and b are uni- 
versal dimensionless constants. Such a linear-in-T form 
is consistent with what has been observed in quantum- 
critical heavy fermion compounds, for both the single- 
particle Green function [3] and order-parameter suscep- 
tibility Q. A linear-in-T relaxation rate signifies that 
the QCP is interacting, i.e., containing a nonzero nonlin- 
ear coupling among the critical modes. By contrast, at a 
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FIG. 3: Scaling of (a) Green's function Q{t,T) and (b) sus- 
ceptibility x{t,T) at the QCP for r = 0.4, To = O.ID and 
Uc{r — 0.4) — 0.085D. The Kondo temperature in this case 
is « 0.029D (for r = 0). For T/D < 5 • 10"^ we ob- 
serve collapse of the data for more than two decades of the 
parameter (7rT)/sin(7rTT), i.e., Qc{t,T) = *(7rr/sin(7rrr)) 
and Xc(T,r) = $(7rr/sin(7rrr)). ^(y 0) oc with 
5 = 0.57(5), and ^{y ^ 0) a t/^~'= with x = 0.68(3). 



Gaussian QCP (whose critical modes do not interact at 
the fixed point), the relaxation rate will be super- linear- 
in-T because the nonlinear coupling itself vanishes as T 
approaches zero 

It is instructive to compare our study with previous 
theoretical treatments of the finite-temperature scaling 
behavior of the pseudogap Anderson/Kondo model. One 
study is perturbative in r, which not only becomes 
unreliable for finite r but also does not allow the study 
of the single-particle Green's function. Another study 
carries out calculations in real frequency at finite tem- 
peratures, but relies on the resummation of a perturba- 
tion series whose validity for the quantum-critical regime 
is not clear Yet another study utilizes a Callan- 

Symanzik approach which requires analytic continuation 
that is problematic as reflected in the noncommutativity 
of the resummation and analytic continuation [l^ ! it will 
be important to check whether that procedure yields a 
Q {oJ,T) that is compatible in analyticity with Q {uj,T). 
As a more specific illustration of our results, we note that 
T^Q (lo/T — ^ 0) is a nonzero constant, which is contrary 
to both the perturbative results of Reference [isj and the 
results of the real-frequency Callan-Symanzik resumma- 
tion for g"{uj,T) 0. 

The scaling of the local correlators in terms of 
ttT/ sin(7rrT) suggests that the boundary critical state 
and the associated boundary operators may be described 
by their counterparts in an effective model with confor- 
mal invariance j3ll |. This is so in spite of the fact that, for 
our problem, the pseudogap form of the density of states 
means that the bulk fermionic component of the Hamilto- 
nian lacks conformal invariance. Hence, our results sug- 
gest an enhanced conformal symmetry that characterizes 
the QCP. 

Summary. We have obtained the full finite- 
temperature scaling functions at the local quantum- 
critical point of the pseudogap Anderson and Kondo 



models. Using the results directly obtained in real fre- 
quency (uj) in the large- A'^ limit, and by showing that the 
imaginary-time local correlators of the physical N — 2 
model have the form of a boundary conformally invari- 
ant fixed point, we succeeded in determining the full 
scaling function in both the quantum coherent and relax- 
ational regimes without using numerically ill-conditioned 
analytical-continuation schemes. We demonstrated that 
the Kondo-breakdown QCP features a linear-in-T relax- 
ation rate for both spin and single-electron dynamics, 
which is consistent with the experimental observations 
in the quantum-critical heavy fermion metals. 
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